clc;
clear;
load x_actual_pcc_robust.txt;
load pcc_robust_xekf.txt;
load pcc_robust_xekf1.txt;
% load pcc_robust_xekf2.txt;
% load pcc_esti_xekf3.txt;

actual = x_actual_pcc_robust(1:60, [1,2,3,18,37,38]);% [1,2,3,18,37,38]
esti = pcc_robust_xekf(1:60, [1,2,3,18,37,38]); % optimal
esti_1 = pcc_robust_xekf1(1:60, [1,2,3,18,37,38]);
% esti_2 = pcc_robust_xekf2(1:60, 1:20);
% esti_3 = pcc_esti_xekf3(:,69:88);

t = 1:2:120;
figure(1)
for i = 1:6
    subplot(2,3,i)
    h1 = plot(t,actual(:,i),'k','Linewidth',1); hold on; 
    h2 = plot(t,esti(:,i),'color',[0/255 161/255 59/255],'Linewidth',1); hold on;
    h3 = plot(t,esti_1(:,i),'color',[253/255,60/255,6/255],'Linewidth',1); hold on;
%     h4 = plot(t,esti_2(:,i),'r','Linewidth',1); hold on;
%     h5 = plot(t,esti_3(:,i),'b','Linewidth',1); hold on;
    legend([h1,h2,h3],'Actual','Robust','Random')
    xlabel('t/min');
    xlim( [ 0, 122 ] );
end

subplot(2,3,1)
ylabel('$C_1^L(CO_2)$',  'interpreter','latex');
subplot(2,3,2)
ylabel('$C_2^L(CO_2)$',  'interpreter','latex');
subplot(2,3,3)
ylabel('$C_3^L(CO_2)$',  'interpreter','latex');
subplot(2,3,4)
ylabel('$T_3^L$',  'interpreter','latex');
subplot(2,3,5)
ylabel('$C_2^G(CO_2)$',  'interpreter','latex');
subplot(2,3,6)
ylabel('$C_3^G(CO_2)$',  'interpreter','latex');

